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ABSTRACT 


CorrelsLtion technique has been used for analysing 
nucleate boiling. A simple experimental set-up has been 
developed to study the boiling phenomena. The analog signal 
obtained from the experimental set up is first converted to 
digital form using analog to digital converter of IBM 1800 
DAGS. The digitized signal is recorded on the magnetic tape 
and the computations are performed on IBM 7044 digital com- 
puter. As a first step, the correlation technique is used 
for obtaining the auto-correlation and povor-spectral density 
for the single gas bubble. Correlation technique is much 
less time consuming than the conventional methods and has 
no probe to disturb the system. This technique can be used 
to investigate the effects of various pa^rameters upon 
frequency response of the nucleate boiling. The qualitative 
effect of bubble multiplicity, the degree of sub-cooling and 
surfaice tension has been presented in this work. In the 
nucleate boiling experiment it is found that the size of the 
bubbles are too small for individual detection and the main 
contribution is duo to the convective currents. The peak 
in the spectral density curve for nucleate boiling occurred 
at low frequency which can be attributed to the convective 
currents. More experiments are needed to make any conclu- 
sions regarding the nucleate boiling process. 



CHAPTER 1 


II'TTRODIJCTIOF 


1.1 INTRO PUCTIOIT 

Development in nuclear reactors and rocket engines 
where exceedingly hi^i heat transfer rates occur in compara- 
tively small areas, have focussed attention on boiling as 
a mode of transferring heat at high heat flux densities. 

To attain these high heat transfer rates by forced convec- 
tion would require excessively high velocities with resulting 
high pressure drop. With nucleate boiling, however, they 
can be reached at much lower bulk velocities. Since there 
are so many variables which enter the boiling process, it 
has been difficult to find a heat transfer correlation that 
will satisfactorily predict boiling heat transfer under a 
variety of condition. Diverse models of the mechanism and 
the means of correlating the heat-transfer data have been 
presented in the literature. Comparisons of some models and 
their associated heat-transfer correlations reveal conflicting 
opinions concerning the mechanisms that explain the highly 
efficient heat transport of nucleate boiling. Pairly detailed 
summarizations of the various nucleate boiling mode2-S are 
presented in references 1 to 5* 

i 

In boiling three different regimes exists nucleate, 
transition and film boiling. The change from one regime to 
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another is a.ccompanied hy marked changes in the hydrod 3 mamic 
and thermal state of the system. These three regimes are 
shown in Fig. 1.1. When the temperatiire of the heating 
surface is helow the fluid saturation temperature heat is 
transferred by convection, forced or natural, depending on 
the system. This region (A3, Fig. 1,1) has been extensively 
investigated and equations have been derived which permit 
the prediction of heat transfer ra-tes. Nucleate boiliag 
(BG, Fig. 1,1) starts when the tempera tirre of the surface 
exceeds the saturation temperature by a few degrees. Adjacent 
to the surface a thin layer of superheated liquid is formed in 
which bubbles nucleate and grow from some preferred spots 
on the surface, such as solid particles or gas adsorbed on the 
surface. Thermal resistance of this superheated liquid film 
is greatly reduced by the agitation produced by the bubbles. 

An increase of the wall temperature is accompanied by a 
large increase of the bubble popula^tion causing in turn a 
sharp increa.S 0 of the heat flux. However, as the temperature 
increases, bubbles become so numerous that their motions 
interact. Under these conditions the nucleate heat flux 
reaches its peak. Transition boiling begins if the tempera- 
ture is increased further. In this region (OD, Fig. 1,1) 
no liquid-solid contact exists. The surface is blanketed 
by an ^^nstable, irregular film of vapour which is in violent 
motion. In transition boiling an increase of temperattire is 
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followed "by a decrease of liea,t flur until a minimum value 
is reached at which point film hoiling starts. This new 
regime is characterised by an orderly discharge of large 
bubbles with a regular frequency and at regulam intervals. 
The mechanism of thermal radiation heat transfer becomes 
more significant as the wall temperature is increased. In 
the film boiling region the heat flux increases with an 
increase of temperature but at a much slower rate than in 
nuclea-te boiling. If attempts arc made to attamn large 
heat fluxes with film boiling, a.s high as those possible 
with nucleate boiling, the temperature of the heated surface 
ms.y become so high as to result in the distraction of the 
heater. This is called a *burnout’. Therefore, the highest 
rates of heat transfer can be s-ttained in the nucleate 
boiling regime, taking care to remain below the 'departure 
from the nucleate boiling* point and avoiding the 'burnout* 
phenomena. Zuber [4] has discussed the regimes of boiling 
and analyzed both theoretically and experimentally, bubble 
growth and the hydrodynamic characteristics of nucleate pool 
boiling. 

1. 2 Review of Nucleate Boiling ; 

Various boiling models have been suggested as the 
reasons for the large heat-transfer coefficients of nucleate 
boiling. The principal models appearing in- the literature 
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is presented below. In each of these models j one heat transfer 
mechanism is considered to be dominant. 

Bubble lAita^tion Mechanism ; 

Considerable experimental eT,H dence shows that an 
appreciable degree of fluid mixing occurs near the heater 
surface during boiling. This mechanism suggests that the 
growth and detachment of the bubbles (Pig. 1.2a) from a 
heated surface induces turbulence that promotes heat transfer 
between the surface and the liquid. The highly turbulent . 
liquid transports heat into the bulk of the liquid. The 
bubble is a passive agent that keeps the fluid agitated but 
is never directly involved in the heat transfer process. 
Schlieren and Shad oW’^T' graph high speed motion pictures were 
made of the thermal boundary layer in the vicinity of the 
bubble by Hsu and Iraham [5]. They found that the development 
of a bubble caused appreciable agitation in the thermal layer. 
As the bubble grew, it appeared to displace the thermal layer, 
and agitation was visible within the adjacent thermal layer 
as far as one bubble diameter from the center of the growing 
bubble. The surface a.rea subject to this agitation is 
defined as the area of influence. 

Yapour-Liquid Exchange Mechanism ; 

The vapour- liquid exchange mechanism (Pig. 1,2b) is 
a means of pumping a slug of hot liquid away from the wall 
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and replacing it with a cooler slug. The cooler slug 
picks up heat from the wall before being transported out 
into the bulk, 1 growing and departing bubble acts as 
the piston in this pumping process. This model resembles 
the bubble agitation model. 

Evaporative Hechanism ; 

A third model (Eig. 1,2c) of nucleate boiling 
places empha^sis on the evaporative process associated with 
the growth of the bubbles. The evaporative mechanism 
accounts for sizeable heat flux rate during nucleate 
boiling. 

The inclusion of the concept of an evaporative mechanism 
at the base of a bubble necessita.tes the presence of a thin 
liquid layer. Such an assumption was made by Moore and 
Mesler [6] in proposing the presence of a micro-la 3 ?"er of 
liquid underneath the' bubble. Experimental verification 
of the presence of micro- layer has been established by 
Hospeti and TIesler [7]. They used a chemical deposition 
technique involving calcium sulphate which contained radio- 
active sulfur 35. The radioactive sulfur made possible the u 
measurement of the amount of deposit residue left by a 
sequential history of many bubbles developing at a particular 
site. Auto-radiographs of the deposits was related to the 
amount of liquid boiled off in generating the bubbles; thus, 



estimates of the thickness of the liquid microlayer under 
the bubble could be made. For average heat fluxes ranging 
from 7000 to 20,000 Btu per hour per square foot, the micro- 
layer thickness was observed to varj^ from 19 x 10 ^ to 

-4 

10 inch. 

Sharp [8] employed three kinds of experimental 
techniques to detect the presence and magnitude of the 
microlayer. The sampler of the techniques was the use of 
an electricaJL continuity probe thait wan inserted into the 
base of a growing bubble. The indication of electrical 
continuity between the probe and the hearting surface was 
interpreted as the presence of a thin liquid layer. The 
more sophisticated method involved the application of the 
optical principles of interference fringes to measure the 
thin microlayer. Sharp was able to make estimates of the 
local thickness of the microlayer a/b various ra-dial positions 
beneath the bubble at several time intervals of bubble growth. 
The meas lured microlayer thickness varied between 4-16 p-in. 
These values are somewhat smaller than those observed by 
Hospeti and Mesler, but it is reraarkable that there is an 
overlapping range of general agreement. Sharp used the 
measured values of the depletion rate of the microlayer to 
predict an average heat flux rate underneath the bubble. The 
average rates varied between 30,000 and 50,000 Btu per hour 
per square foot for the area under the bubble. 
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Voutsinos [9] iias investigated tlie growth and 
evaporation of the microlayer underlying a huhhle forming 
on a glass heater surface using laser interferometry and 
high speed photography. The results for the single bubble 
indiceted thsf the microlayer thickness is of the order of 
5 P”in. Subsequent analysis of these results confirms that 
the microlayer evapora.tion phenomena is a significant heat 
transfer mechanism, representing approximately 25 percent 
of the total nucleate boiling heat transfer rate. 

Grraham and Hendricks flOj have proposed a model of 
nucleate boiling where time and surface area averages of 
the following basic heat transfer mechanisms have been taken 
into consideration (i) transient thermal conduction through 
the liquid thermal layer in the vicinity of a nucleation 
site that is preparing to bear a bubble (ii) evaporation 
from a microlayer surface underneath a bubble that is attached 
to the heater surface (iii) turbulent-free con-vection that is 
taking place over the surface areas not covered by bubbles 
(A zone of enhanced convection occurs in the vicinity of 
a growing bubble). The overall model presented by Graham 
and Hendricks emphasizes that nucleate boiling should be 
thought of primarily as a series of transient heat-transfer 
processes over distinct regions of the surface, 

Mikic and Rohsenow [11] assumed that the main 
mechanism of heat transfer in nucleate boiling is transient 
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her-t conduction and subsequent repls.cement ■ of the superlieated 
layer around boiling sites associated with the bubble depar- 
ture. They found that the lioat flux versus wall superheat 
relation depends on the cavity size distribution of the 
boiling surface. Torikai anc Ya,majzaki [12] studied boiling 
heat transfer mechanism by using the electro-conductive 
glass plate and high spoed motion picture camera. They found 
throe mechanism of boiling heat transfers First mechanism is a heat 
conduction through the thin liquid film in contact area 

of bubbles on the heating surface, which is main part in 
heat transfer at the high superheating temperature in a 
wotta.ble surface. Second mechanism is a heat transfer in 
natural convection by bubbles and in forced circulation of 
flow, which is dominant a.t low superheating temperat’ero in 
any surfaces in saturation of liquid. Third mechanism is a. 
heat conduction in unsteady state by contact of liquid on the 
hea.ting surface during the detachment of bubbles, which is 
large part in heat transfer in sub-cooling liquid. Tong [13] 
has described the boiling raechanisms at various flow conditions 
and has given emphasis to the- prediction of incipient nucleate 
boiling, effects of surfaC'v. fouling, and transient cooling. 

G-anic and ifgan [14] studied the bubble growth jn 
pool boiling by following the growth of the bubble by the 
use of high-speed photography and simultaneously measured 
the temperature in the bubble and its liquid environment by 
a fast response thermocouple. lUhey measured the continuous 
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temperature distribution around a single bubble by changing 
the position of the thermocouple iprobe and analyzed a large 
number of bubbles. 

Fail [15] used a scanning Electron Microscope (SEM) 
to examine pool boiling nucleation sites on 504SS in contact 
with distilled degassed water, ilucleation sites were 
located during the boiling process 13?- noting the coordinates 
of the bubble locations on the surface. He mB.gnified these 
locations using the SEM and made photographs of nucleation 
sites. The radii of the cavities vanied from I4.8 p-in to 
180 |i"in, and width of groves varied from 75 ii-in. to 
154 (i-in. 

In experimental program of transient boiling under 
low pressure condition was initiated by Schultz et al [16] 
to study the single artificial cavities of known geometry 
and size. The technique employed to characterize the tempe- 
rature field relies on experimental heating period of such 
short duration that convection effects can be neglected and 
the temperature distribution computed analytically using a 
pure conduction model. They examined the role of the develop- 
ing transient thermal boundary layer on the stability of the 
just nucleated vapour bubble throu;^ the use of the analysis 
of hi^ speed close-up films of the growing nucleus. 
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Rogers and Mesler [17] used a special surface thermo- 
couple and measured large temperature fluctuations during 
nucleate boiling of water on chromel P surface. They 
found that the temperature dropped 20 to 30°P in about 
2 msec, but required 10 to 20 msec to recover. They 
postulated that surface is cooled during initial bubble 
growth by evaporation of a microlfyex into the bubble. To 
test this hypothesis they developed a technique .to photograph 
a bubble growing from an artificial site as the surface temp, 
was measured. They found that the result is in consistent 
with their hypothesis. 

Theofanous and Patel [18] have observed that the 
initial to final vapour density ratio (when significantly 
greater than one ) strongljr influences vapour bubble growth 
rate in the inertia and the inertia/heat transfer regimes. 

Since boiling heat transfer is controlled by generated 
bubbles on the heating surface, to understand the nucleate 
boiling a better knowledge of bubble is desirable. Many 
aspects of the study of bubbles, such as formation of 
bubbles, bubble shape, growth rat e, motion of the bubble 
maximum size, and forces acting on a bubble have been studied 
by many investigator [19-25]. By far the most frequently 
employed technique in the study of boiling has been the 
analysis of hi^ speed motion pictures, A frame by-frame 
analysis of the pictures provides the two-dimensional 
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projected size and the loca,l velocity of the buhhles. 

The main disad vantages of this approach are that the 
analysis is very time-consuming and yields only a tT:'XO- 
dimensional representation of the hu'b'ble. 

1.3 Gorrelation Technique ; 

Correlation procedure takes tvo functions of time x(t) 
and y(t )5 delays one of the functions by an amount x, then 
multiplies them together and averages the result over some 
finite time. If x(t) = y(t) an 'autocorrelation* fTHiction 
is involved j if x(t) and y(t) are different, a * cross - 
correlation', function is obtained. When x(t) =; y(t), the 
resulting auto-correlation function (t) is symmetric 
about T = 0, with O^^(O) greater than or equal to any other 
value of In practice, thez-e is usually a sharp 

maximum at O^^(O) so that one is able to distinguish G^f(O) 
clearly from any other value. In accordance with the 
Wiener theorem, the power daisity spectrum is obtainable 
as a cosine Fourier transform of auto-correlation function. 
This is the only means by which power density spectrum of 
random function is cal-culated. Thus, the correlation 
technique is the application of correlation to the spectral 
analysis of the random phenomena. This technique can also 
be used for the spectral analysis of periodic and aperiodic 
function (Ref, 26,27). 
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Sines nucleate boiling is a random process, tiie 
techniques of noise analysis should be usefiil in the 
study of the process. The first iuTestigation of boiling 
water nuclea.r reactors utilizing noise analysis techniques 
was performed by irgenne National laboratory [28]. Transfer 
functions, spectral moments and correlation functions were 
obtained for the BORAS-IV reactor. 

Mcilnight and Ra,m [29] used correlation technique 
to determine the power spectra and auto-correlation for a 
single gas bubble rising in water and for a nucleate 
boiling system. Both analog and digital method of getting 
power spectral density were used. In analog method, the 
signal was squared and integrated over a measured time 
intervo.1 and this information was used to calculate the 
power spectral density. In digital method, the auto- 
correlation curve obtained by using a signal correlator, 
was digitized and the Fourier analysis was performed 
using these digitized points. These functions were then 
analyzed to determine the bubble size, velocity and 
frequency of generation. They also showed that how the 
hes-t fluK, the degree of subcooling and surface tension of 
the liquid affect the power spectrum of the nucleate boiling 
system. 

McOurdy [30] showed the relationship between the 
spectral density of the spatially -averaged temperature 
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(of a heated filament) and the auto-e^nd cross-spectral 
densities of the temperature induced at individua.1 sites. 

By analyzing the trace of the spa,tially -averaged temperatures 
of the system, he obtained its spectral density. Using 
correlation technique , he investigated nucleate pool 
boiling in saturated and subcooled ws-ter in normal and 
low gravity fields and obtained the power spectra of the 
temperature fluctuations at the nucleating surface. The 
spectra of the surface temperature fluctuations were 
found to contain a number of peaks, occuring- at the same 
frequencies vrith which bubbles were emitted from the 
various active sites on the heating surfooe, as determined 
from the hi^i-speed films. The spectrum of the temperature 
fluctuations was markedly influenced by the degree of 
subcooling and by an inverse gravity field. McCurdy found 
that generalized harmonic analysis is a useful tool for 
the study of heat t ransfer phenomena which possess a random 
character such as the temperature fluctuations which occur 
in nuclea.te boiling. 

Osermely and Ostrander [31] have used the correlation 
technique for determining the frequency response function of 
a heat exchanger. They used a concentric tube water-to- 
water counter flow heat exchanger which was subjected to 
temperatuia perturbations at the tube-side inlet. The 
perturbations were produced by a discrete interval binary 
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noise signal# Correlation and spectral analysis for 
obtaining tbe frequency response was accompanied by means 
of IBM 1620 digital computer with the tube inlet and shell 
outlet temperature histories on the punched paper tape as 
the computer input. The frequency response function 
determined by the correlation technique was found to be 
in reasonable agreement with the frequency response function 
obtained by the more conventional step and pulse-input 
technique. 

’Determination of the frequency-response function 
of a system by correlation technique has the advantage 
over the other techniques that interference of signals with 
noise is reduced in the end results. This technique should 
be used whenever noise is a problem. Correlation technique[ 29] 
can be used to obtain infoimation previously obtained by 
using time-consuming motion picture methods, specifically, 
that the two-dimensional projected size of the bubble, its 
local velocity, and the size and frequency distributions 
of the bubbles can be determined, 

1,4 Present Work ; 

The present work consists of utilizing correlation 
(Noise) techniques for detecting and analyzing nucleate 
boiling. Correlation techniques have the advantages of 
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i) Consuming less time for analysis, 

ii) Fo probe to disturb the process, 

iii) Easy measurement of tbe power spectrum 

iv) Different paremeters affecting the process can be 
correlated . 

v) High degree of precision even in the face of severe 
input noise to signal ratio. 

Correlation techniques thou^ bCvSc-d on simple 
principles, are usually complex in the sense that digital 
computing machines or special electronic circuits are used. 

Here, as a first stop single ga.s bubbles are 
analyzed for their auto-correlation and power spectral 
density and later nucleate boiling is studied. 

In the present work, correlation techniquo uses 
digital method for calculating the auto-correlation and 
power density spectrum of the signal. The digitization 
of analog signal has been accomplished by an anaHog-to- 
digital converter (ADC) using IBM 1800 DACS. The digitized 
data has been recorded on the magnetic tape. The compu- 
tations have been carried out on an IBM 7044 digital 
computer. Output from IBM 7044 are recorded on magnetic 
tape, and this tape has been later brought back to IBM 
1800 for final plotting of the results. 
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Chapter 2 of the thesis deals with the hackground 
theory of the correlation. Here the computation of the 
correlation fimction from experimental data on the "basis 
of time averaging has been discussed.. Also described in 
this Cha,pter are the alterna.tive methods for data acquisi- 
tion and analyzing the analog signal. 

Chapter 3 contains the description of the physical 
layout of the system, system hardware and software. In 
the ■ software , the functions of the varioua -computer 
program and sub-routines used for recording and analyzing 
the data have been discussed. 

Chapter 4 describes some of the results obtained 
by this method for single bubble as well as some prelimi- 
nary results of nucleate boiling. Results of the surface 
tension variation on single bubble as well as nucleate 
boiling have been presented. 

Chapter 5 concludes with some recommendation for 
future work, especially to study the nucleate boiling 
phenomenon. 



CHAPTER 2 


THEORY 

2 . 1 Discussion of Analj/tical Method ;, 

In manj/ areas of science and engineering, the 
accuracy of experimental measurements is limited by the 
signal to noise ratio of the system. This fact has led 
to the development of the field of information and statis- 
tical communication theory. Harmonic analysis technique 
has been generalized to include randomly varying signals 
(noise) with the aid of Fourier transform. This generaliza- 
tion is quite important because random signal or noise, 
contain information describing the character of the system 
from which they emanate. Random signal analysis (noise 
analysis) has been used in many diveirse fields, such as 
nuclear engineering, communication, radio astronomy, 
seismology, neurology etc. The brief treatment of random 
signals (noise ana,lysis) which is included in this chapter 
follows essentially that of Thie [32]. 

Correlation Functions and Power Spectral Density Functions ; 

Consider some physical process which gives rise to 
a time-vaiying signal f(t). The signal may be simple or 
complex periodic or have the character of noise, that is, 
be random varying. Auto -correlation of a signal f(t) is 
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o'btalned by telcir-£; til© tirae aver?, 3 3 of f(t) fft-'x), Fbero 
f(t"T) ig came signal d3la3''GC by time t. Matliematically 
the auto-correlation function -r ) is defined as ^ 

T 

t= lira / f(t) . f(t~i:) dt (2,1) 

where I j_,3 period of the signal. 

Important properties of suto-correlption function 

f ( ) : 

1 ) 


2 ) 


is sn even fimction of t i.e,. 

hen maximum value for zero delajr time (t = O). 

i.e, C^.p(0) ^ ^f f ^ 

is independent of the time origin. Ihis 
means that auto-correlation function of f(t) is the sane 
s-s ths.taf.: t-'t^) where t^ may hanre any value. 

]?or rsal f(t)’S; O^^('r) is a real function of t. 

Applying property ( 1 ) of the a-uto-correlatioio 
xuncxion, equation (2.1) can be written as, 


T 


^^^(t) = lim /■ f(t).f(t+i:) dt 

ip ^ 00 _p 


( 2 . 2 ) 


for the case in which the signal f(t) is given as f discrete 
data points (electronically by digitizing an analog signal, 
the auto-correla^tion function may be determined by 
replacing the integral in eqnation (2.,2) by a summation as, 
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C 



wr y ^ ( "fc -j ) 


(2.3) 


The normalised auto-correlation function is defined 
with respect to deviation from the mean and is given hy, 


(/f^( r) 


Lim — i — J ff ( t)-f ] . [f ( t+T)-f ] dt 

T->“ 2ai -T &T 

X 


lim —7 ^ — 
N->-" H 


N 


y [f(t .)-f] .[f(t .+T)-f] 
j=l ^ ‘J 


(2.4) 


where f is the mean value of the time -varying signal over 
one period of a periodic signal or over the discrete 

2 

number of digitized values of a ra.ndom signal, and cr^ 
is the variance of the signal f(t). Strictly speaking, the 
auto-correlation function is compiited in the limit of 
N 00 or T “ , but in practice, however, one treats a 

finite amount of data and computes 0^^ for r’s upto 

T = MA t < H At 
m 


for S' data points spaced by At. 

IT-t/ A t _ P 

(2.5) 

whe re , 
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md . 


1 f 

^ jii ' 


I '-‘j+C-r/tt) ' (‘-6) 


N-M 


3=1 


The auto -correlation function G^^(t) and the nomalized 
auto-correlation function are related by, 

C^^(t) = ct| + (f)^ (2.7) 


Similarly, if two signals, f(t) and g(t), either 
an input and an output or two outputs to a known input, are 
obtained from a system, and one is delayed and the average 
value of the product of signals is determined, then the 
result is the cross-correlation function of the two signals. 
This may be written as, 

1 ^ 

Of (t) = lim / f(t).g(t+T;) dt (2.8) 


As with the auto-correlation function, one can define the 
normalized cross-correlation function as, 

T 

(t) = Lim / [f ( t )-f ] . [g(t+T)-g] dt 

^ T-><» 2o. 0 T ~T 

^ S (2.9) 

and for N synchronized data points spaced At apart 

for both signals, the normalized cross-correlation function 

may be calculated for t upto = MAt < lAt as, 
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{2^^g(T) = [ 


1 




N-T/At In 
J ^ 


j +( xf LX) 


- nihf «g 
(2.10) 


■whe re 


-1 

5 = I :L 


o 1 ^ — O 

4 = I .1^ (f,-f)^ 


The cross— CO rrelatior function can he described as 
representing the degree of confornitj'" between two signg^ls. 

The correlation functions are useful in describing 
a system' s response in the time domain. The impulse 
response function h(t) of a system is related to the input 
f^( t) and output fQ(t) by the following convolution integral, 

f (t) = f f.(v) h(t-v) dv ( 2..J.1) 

o J 1 

..00 

Using this relation, the correlation functions defined in 
equations (2.2) and (2,8) can be wn-itten as, 

O^^(t) = / h(v) C^^Cr-v) dv (2.12) 

«-00 

where Cj^j^(v) and ^j_Q('i^) auto— correlation of the input 
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and the cross “Correlation between the input and the output 
respectively. If the input is purGl3r random^ then 
G..(t-v) is a delta fimction, xfhich indicates tha,t the 
cross-correlation is relented to the impulse response 
function h(v). 


¥ ic ner The o reni for Auto-Go rrelat ion s 


Wiener theorem for auto-correla,tion states that auto- 
correlation function of a random function 3Jid the power 
density spectrum of the random function are related to 
each other by a Fourier cosine transformation as given by, 


C( t) 


where P(f) 
P(f) 


f P(f) exp (iw-r) df ^ 



P(f) cos coT df 


O 

lim ^ lg(f)l^ 


Lim 

/ 

C(t) 

exp 

^ ^ 00 
m 




Lim 

C\J 

G(t) 

cos 

T -^00 

m 

0 




( -iw-r) 


wT dT 


o = 27Cf 


dT 


(2.13) 


(2.14) 


With W discrete values of G(t) spaced At apart, the power 
density spectrum can be expressed as(Ref* p9‘)j; 


P(f) 


ff + I 


M 

y C!(jAt) cos 2% f ,]At 

j=l 


M 

= At C(o) + 2At y G(jAt) cos 27 i:fjAt (2.15) 

j=l 
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In the computer program described in the next 
chapter, eqn. (2.15) has been used for calculating the 
power density spectrum. 

2.2 Data Recording and Analysis ^ 

There are two general method', of recording and analysis, 
of the data. 

1) Ana,log or continuous and 

2) Digital or Discrete 

'The former may involve the use of analog computers 
if done off-line. If done on-line, the term continuous 
is appropria^te because the electronic suialyzers use the 
filters, amplifiers, and integrs.tors that handle conti- 
nuous function of time. Digital da,t a handling is invariably 
done off-line because a. digital computer is involved. It is 
quite possible that analog and digital methods may be used in 
a single experiment. 

Continuous recording may be accomplished by standard 
chart recorders or by magnetic tape recording. Having once 
recorded an experimental data on ma.gnetic tape, it can be 
rerun and reansHyzed at any time and as often as desired. 
Data so recorded is easily stored and is compatible with 
both analog and digital methods of analysis. Also by using 
different recording and playback speeds, one can achieve 
various time transformations in the analysis. Recording 
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continuous signals on magnetic tf.pe may le done in two ways 
a) Direct amplitude recording in ■'wliicli the degree of 
magnetization is proportiona.1 to the signal , (h) .Frequency 
modulation of a carrier directly recorded at constant 
amplitude. Ihough the first has been highly successful 
in audio frequency ( 20 cps to 20,000 cps) application, the 
second is used in data recording oGca,use frequency 
response to low frequency is needed. 

When digital analysis of a continuous signal, is 

desired, an analog to digital conversion must be perfoimed. 

The rate at which the conversion of analog signal to 

digital signal is done is decided by the sampling theorem. 

The theorem states that if a low pa.ss band limited signal 

of bandwidth fjjj- is sampled at the rate greater than 2fj;j. 

samples per second, then the original signal can be comple™ 
determined 

tely/from the sampled signal. Sampling time At is then 
related to the limiting frequency as, 



or, 

f _ -1- = 1 ^ X 

2At 2( digitizing interval) 2T 

i 

where the sampling frequency % is called the ITyquist 

i' 

frequency or cut-off frequency. Frequencies above f-j^^j cannot ; 
be detected and therefore the numerical integrations have 
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as th.e highest limit. In the event frequencies above the 
Hyquist freqv.ency sare actually present in the continuous 
signai-, then a phenomena, known as ’aliasing’ occu.rs, that 
is, in the digital analysis these higher freauencies will 
appear below the Hyquist frequenc“" and will be indistin- 
guistable from the lower frequencies. To' avoid difficulties 
from this source , ono should choose f so that it exceeds the 

o 

maximum frequency for which the spectral density is signi- 
ficant, If fj^j is this maximum frequency],' then requiring, 


's i % 


means that the ssanpling interval 


^t should satisfy, 


At 



1 


2f 


M 


In words, the time interval between successive samples 
should be such that the sampled da,ta contain at least two 
samples per cycle of the highest frequency. In practice, 
ten to twenty/ samples per cycle of this highest frequency 
are used when noise is present. 

If the auto -correlation function is determined upto 
= MAt < HAt, then the higgles t frequency will be 
fjvj = E/2z.^. This means the spacing of the frequencies 

Af = 1/2t,^ 01’ 1/2M At . However, the spectra,! resolution 

in 

due to a finite duration of sampling, T, is equal to 1/T 
or This limitation results since the auto -correlation 
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function is not known for x > Because of this one 

defines a spectral lag window for infinite t’Sc 


h(T) 1 1 

for 

!t 1 < X 

' ' — m 

= 0 

for 

l'>:| > •'m 


There are a variety of spectral windows, h(T), which 
can he used to minimize the effect upon P(f) of contribution 
outside the range f + Af. 

Two frequently used windows are, 


Kanning 

lieroming 

Since the power spectra obtained by eqnation (2.15) are 

spaced by 1/2 Tj^, one of the above windows is used for 

smoothing them. In present work, power spectral density 

and the smooth power spectral density are calculated by 

using the following expressions [29], 

M 

P(f) = At C(0) + 2At y G(jAt) cos 2Ttf .At (2.16) 

j=l 


0.5 + 0.5 cos 


nx 

''-m 


0.54 + 0.46 cos 


%x 

'^m 


P' (f ) = 


0.25 P(f - 



+ 0.50 P(f) + 0.25P (f + 



) 


(2.17) 


Equation (2.17) has been used to write computer program which 
calculates smooth power spectral density from power spectral 
density obtained by the use of equation (2.15). 



CHAPTER 5 


SXPEPIMEHTAI SET-ITP /JTP jiITALYSIS 


3 . 1 Sy stern Peso rip t ion g 

The complete inyestlgation procedure cards broadly 
be divided into two main parts: 

(1) Hardware Development , and 

(2) Software Development 

Before describing details of the experimental set-up 
and data- logging procedure, a blocl: diagram representation 
of the system will be outlined. S; 7 stem representation is 
given in Pig. 3- » 1' whe re all the functional elements of the 
complete scheme are depicted. 

In general, the p 2 X)cess under investigation does not 
yield information which can be processed directly without the 
use of transducers. The transduCvir converts the message to 
an analog signal, a time varying olectrical quantity, such 
as voltage or current, which is better suited for further 
processing by the data processing system. 

For convenience and clear understanding, these things 
have been isolated as distinct identities, though in an 
actual system the separation may not be so obvious. Also 
shown are some of unwanted noise or disturbances which are 
inherent in any set-up. 



Thougiis the description sjic; details of each iKn 
block will he gi-ven separately j a look at the block diagram 
gives the signal flow information, 

A parallel beem of white light (Block 1 , Bi;. 
passes through the transpa^rent gla„S3 of the boiling 
(Block! 2, ^ig. 3.1), having a bent copper tube with a 
hole in it to produce the bubble by passing g8.s or bubbles 
generated from a. heated surface. The light beam is attenuc!,ted 
by bubbles rising through the medium, fluctuation in the 
light signal caused by scattering of light by bubble is 
detected by a photo:^ultiplier tube which acts as an input 
transducer (Block 3, Fig. 3.1), converting the light signals 
to electrical pulses. Linear amplifier (Slock 4, Fig. 3*1), 
amplifies this weak signal and transmits through coaxial 
o-able acting as the transmission medium (Block 5, Fig. 3.1), 
to the interface equipment of the data processor. Trans- 
mission mediura also sometimes called channel is an, electrical 
connection in between the linear amplifier and input interface 
(Block 6, Fig. 5.1) bridging the distance between the source 
and data processor while maintaining the true characteristics 
of the signal. In this case the transmission medium is a 
coaxial cable of 75 ohms total impedence. Transmitted signal 
is fed to the input interface and filter where unwanted 
noise signals picked up during transmission are suppressed 
and the output of the filter is connected to solid state 



multiplexer of IBM 1800 (Block 7, Big. 3.1). This analog 
signal is converted into digital signals "by an analog to 
digital converter (IBG) within the IE!! 1800 system. 

Digitised signal values are recorded on the magnetic tape 
of IBM 1800, using softxfare to he discussed later. 

Selective samples of data are processed on IBM 1800 for a 
quick analysis of the data and the majority of the samples 
have "been processed by IBM 7044 (Block 8, Big. 3.1). Signal 
was analysed for its auto-coirelation and power spectral 
density. Since there was no plotter available in IBM 7044, 
output was put on a magnetic tape and brought back to IBM 1300 
for final plotting of the results. Detanls of the hardware 
componenas is described below: 

3.2 Hardware Development : 

The experimental set-up consists of the following 
components j 

(1) light source 

(2) Vfater tank 

(3) TMT and linear Mplifier 

(4) Transmission Medium 

(5) Input Interface and Filter 

light Source ; 


The optical set-up details has been shown in Fig. 3. 2. 
A 12 V, 3 mA tungsten filament bulb was used as the light 
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source. Power was supplied by a 30 volt power supply 
unit which was adjusted to give the desired voltage. 

A reflector was provided to focis the beam in one 
direction and avoid scattering of light. Phis light 
beam was passed through a G-.I. pipe of diameter 4.2 cm 
and at the other end of the pipe a convex lens of focal 
length 25.7 cm was placed, lis'tance between the light 
source and lens was so adjusted that the light source 
falls exactly at the focal length of the lens to give a 
parallel beam of light. The who-le set-up is enclosed 
in an aluminiimi box covered from top to protect from 
extraneous 3_ight. The beam is collimated through a 
hole of 6.5 mm diameter in the front to allow a narrow 
beam of light comparable to bubble size to come out of 
the box. 

Water Tankj 

Since experiment is perfor'^ed in two parts (a) with 
gas bubbles (Pig. 5.2a), (b) nucleate boiling bubbles 
(Pig. 5. 5), two different test containers are used. 

Test Container for G-as Bubbles 

Pig. 5.5 shows the pictorial view of the container. 
This consists of a simple ^ass jar of 57.5 em x 8 cm. This 
jar is filled with demineralised water to a height of 26 cm. 

generated in the jar from an orifice in a 


Bubbles are 
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copper tube submerged in the iar. An air cylinder was 
used to pass air through the copper tube. Bubble multi- 
plicity and frequency of emission is changed by varying 
the flow of gas through the orifice. Glass container is 
wrapped in black paper to protect from unwanted light 
except for small portion to transmit the light signal. 

Water Tank for Nucleate Boiling; 

In the study of nucleate boiling the bubbles are 
formed by hearting a copper rod pointed at the tip. The 
test piece is pointed to have preferably single nucleation 
site. However, this goal- is partially achieved during the 
experiment, since it is observed that nucleation is 
occurring on the surface other than the tip. Fig. 3.4 
shows the dimensions of the test piece. This test piece 
is heated externally by nichrome heater ribbon. To have a 
good thermal insulation, the copper piece is, passed through 
a push fit teflon plug before fitting it at the centre of the 
water tank. Heat flux variation is achieved by varying the 
voltage across the nicrom ribbon. 

-Fig, 3.5 shows the constructional details of the 
water tanJi. The water tank filled with demineralized 
water to a height of 11.5, cm is provided with glass windows 
on two sides for allowing the light to pass through the tank. 
The dimensions of the water tank are 45.5 cm x 30 cm x 44*5 cm 
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whicli is very large compared to the 6.35 mm dia. of the test 
piece to minimize the effect of convection cxirrents. Heating 
of water in the tank is achieved, by two external heaters of 
1 Elf each, mounted at the bottom of the tank on either sides 
of the test piece. Different constant bath temperature is 
achieved by varying the supply voltage to the heaters by 
means of a variac. It has been observed that for a given 
water level in the tank, and a given voltage to heaters, a 
particular constant maximum bath temperature is obtained. 

A series of experiments with vany'ing heater voltages are 
performed keeping the level of water constant, and when the 
temperature rea.ched a steady sta:.te value, it is noted down. 
Voltage V3. temperature graph is plotted and this information 
is used in conducting the nucleate boiling experiment. The 
temperatures are recorded using a alcohol thermometers, the 
bulb of the thermometer is kept at 1.5 cm away from the 
tip and very near to the surface of the tank. 

Photomultii3lier Tube ? 

Photomultiplier tube (PMT), is the most important 
component of the whole system. This senses the fluctuation in 
the li^t signal and converts it to electrical signal. PMT 
is essentially a cuirent source. 

PMT used in this experimental set is RCA 68101 . This 
is a 14-stage, 2" dia. head-on type of photomultiplier tube. 
The photocathode material is 0s~Sb., 6810 a features hi^ 
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q_u8-iitum efficiency 5 fast time resolution characteristic, 
low dark current, and high current amplification charac- 
teristic. The spectral response of 68101 covers the 
range from about 3000 to 7000 angstroms (i). Maximum 
response occurs at approximately- 4000 i. Therefore, it 
has high sensitivity in the blue and less sensitivity in 
the red regions of the visible spectrum. 68101 is 
rela-tively insensitive to the effects of extraneous magnetic 
and electrosts,tic field. Maximum operating voltage for 
this tube is 2400 volts DC only. However, in our case it 
is much below this maximum rating. Voltage divider arrange- 
ment and basing diagram has beoxU shoxirn in Fig. 3.6 and 
3.7. 

linear -Implifier; 

Since the output of PMT is very small, it become very 
essential' to amplify the signal before coupling it onto the 
transmission line, 

Implifier described in Fig. 3.8 consists of -3 stages s 
(i) Operational amplifier (a) is input buffer having unity 
gain. Input impedence of this stage is kept very high to 
match the impedence of the PMT, 1 gain of 1000 is achieved 
by 2 stage amplification (b) and (c). Transistor and Q 2 
are current amplifier to drive the coaxial line having very 
low input impedence. The overall band-width of amplifier 
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sij3.igs is kep'G upto 1 S.Hz as the signal expected is in 
low frequency region. After amplification, signal amplitude 
is in the range of 0.25 - 3.5 volts to ho compatible with 
ADC input of IBM 1800. 

Transmission_ Me_d _i-umj. 5 

The outpi.it of linear amplifier is transmitted to 
IBM 1800 by means of a coaocial cable of 75 ohms total 
impedence. Coaxial cable is chosen to minimize the noise 
pick up during the transmission, since the signal carrying 
cable is quite often run close to the source of main voltage 

Input Inte rf a ce and Filter : 

In course of elec trie a, 1 sigrxil transmission certain 
unwanted and undesirable effects take place, namely attenua- 
tion, distortion, interference, and noise. Attenuation 
reduces the signal strength. In tiiis case attenuation 
is very negligible, since signal ''is transmitted over a 
very small distance. Distortion is signal alteration due 
to imperfect response of the system to the desired signal 

itself. Interference is contamination by extraneous sig-' 
Foise is 

nal. /random and unpredictable electric signals which come 
from natural causes both internal and external to "the system 
when such random variations are edded to an information 
bearing signal the information may be partially masked or 
totally obliterated. 
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Keeping above facts in viei/^ Input interface a,nd 
filter is designed to suit tbe signal at tiie receiving end. 

It consisted of txfo stages as shown in Fig. 3,9. 

First. Stage; acts as input interface. High input imped ence 
unity gain amplifier is used. Matching is achieved by 
resistance network at input of amplifier. 

Second Stage ; acts as a tunable filter, to cut off all 
signals of 50 c/s pick up and other undesirable frequency. 
Approximate values of resistors a.nd capacitors are indi- 
cated in Fig. 3.9. 

A short description of IBM 1000 is given before giving 
the details of data logging progra.mme. 

I BM 1800 Bata Acq uisition and Control Simulation (DAOS) ; 

IBM 1800 DAOS is capable of accepting one or more 
analog inputs. Analog signals representing the physical 
phenomenon are to be amplified and fed to the analog to 
digital converter (ABO) of the ISOO.. When more than one 
signal is to be processed simultaneously a multiplexer 
scans the signals sequentially. There can be as many as 
10 24 analog input teiminals attached to one ABC. Out of 
these, only 256 terminals can be solid state. 

There are two types of multiplexer available. One is 
Relay multiplexer and the other is Solid state. The Relay 



45 


niul'tiplsxGr c-ccepts an analog signal wliich. is dif isrsniial 
in nature and is within + 5 volts dynamic range. The 
multiplexer scanning frequency from one signal to another 
is 100 cps maximum. Solid state multiplexer can he used 
for high speed data a.cquisition. The signal attached to 
these terminals should have one terminal grounded, and 
accepts signails of + 5 volts. Solid state multiplexers 
have a scanning frequency of 20,000 points per seconds. 

The ADC of IBM 1800 has the ahility to converting 
bipolar analog signals within + 5 volts range to digital 
va,lues . It includes a buffer amplifier and has program 
selectable resolutions of 8, 11 and 14 bits. The larger 
the bit value the hi/^er the accuracy and more conversion 
time. The conversion time for 8 bits, 11 bits and 14 bits 
resolutions ane 29 p sec, 56 p soc, and 44 \x sec. respectively. 

The input impedence of ADC is of the order of 
10 megaohms. 

Data word of the i^DC; 

The data word developed in the ADC is compatible with 
IBM 1800 word format. Each word starts with a hit alloted 
for the positivo/negativc sign an<2 ends with a bit reserved 
for the overflow indication. The output of the ADC is in 
the binary form which is stored in the memory. 
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TSX 8.nd MPX Exeauti-ge Profi:rams; 

Por Ij-jM 1800 DAOS, tliere are two programming systems 
available viz.. Multiplexing iSsocutive (MPX) and Time- 
sharing Executive (TSX). TSX program divides the entire 
core area into -tvro parts, the skelton and V-core. In the 
Skelton area, all the importa.nt system programs and sub- 
routines and a few frequently used user's program axo 
permanently stored. The skelton area occupies approxi- 
mately 8 to 10 thousand words of memory. IBM 1800 has 
16 K memory. Therefore remaining memory (about 6 K) is 
availa,ble to the user under TSX system. 

The MIX program divides the entire memory into 
several segments knoxm as Areas. All the important system 
prograjiB and subroutines are permanently stored in the 
skelton areai, occupying approximately 12 K words. Best of 
the memory is divided into segments to execute programs of 
different execution priorities. MIX requires large core 
memory, at least 24 K word. The MIX system, generated at this 
installation with bare minimum facilities takes 11100 words 
of core, leaving only around 6000 words for users. 
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3*3 Software Development i 
Data Loggings 

The data logging program is given in Appendix A. 

The data logging progra.m consists of suh-progra-ms which 
can he modified to give different sampling rates and variable 
number of samples. The programs are also flexible and 
new subroutines ca,n be added and/or deleted, if deaifed, 
to make the da.ta logging program very versatile. 

Since the maximum frequency components of the 
bubble signal expected is around 40 - c/s, a sampling 
rate of 80 would ha.vo served the purpose. But the sampling 
rate chosen for data, logging is <50 sample /sec. in this 
case. Signa.1 is sampled for 8.152 sec. to give 2048 samples. 

The digitized data in the binary form is finally stored 
on a magnetic ta.pc. . Since there is only one tape unit 
available some opera.ter intervention becomes inevitable 
for changing; the magnetic tape. 

The data logging program is a link program. The 
link is needed to accommodate the main program and- sub-programs 
in the memory of IBM 1800. 

There are two links for data logging program. 



FIRST LII'TR ; ETIGI 


j/Tl'G-l (ma,in program) 

Si'JTP 

ZSRCH 
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ivi.'lliji! 

TSRl/F- 

I'/RIT 

SSCOKP LIIK OUTPT OUTPT (main progrejn) 

IPIOG 

I/IIT 

DTLG-1 i The main progrcni of the first link controls the 
sequence of operations to he performed. I)€;pending on the 
status of 16 data switcLus (PATSW), control is passed from 
one subroutine to another. 

SAMP: This subroutine is most important among all sub- 

routine because this decide the sampling time and starts the 
A/D converter to acquire data of a specified record length. 

Z3RQH : This subroutine- is executed only once in the beginning 

of the programme. This subroutine searches for the last logged 
data on the tape and makos the tape ready for recording 
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more data. 

data, logging program Iceeps track of time, date 
and month. This is achieved by the subroutine KAIBN. 

Time, da.te and month must be entered when executing the 
data logging program. During the logging period the 
subroutine ICALEh keeps .updating the time, date and month. 

SUI'-IRY ; At the end of data logging, a summary of what has 
been put on the tape is veiy ranch essential. SDl'IRY 
subroutine is cailed after the end of the data logging, 
and this gives the serial no. of the data logged, status 
of the data switches and the time. Data switches represent 
the information regarding the condition of the system, 
number of data points, gain of the system etc. which 
is recorded on the tape. 

TSA'VE i This subroutine serves the function of unloading the 
tape to save the recorded information. 

WAIT ; This subroutine comes in action whenever the data 
switch 3 is put in UP position. The data logging is 
suspended till the data switch 3 is in UP position. Since the 
UP position of data switch 3 can keep the computer in wait 
position, the experimenter can chenge the signal condition 
at the transmission end, without stopping the data logging 
program. As soon as the data switch 3 is put DOW, the 
data logging once again starts. 
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OUTPT ; This main program of the second link, is called 
when plotting of the signal is desired. A typical signal 
from the tape can be selected for plotting. OUTPT selects the 
roquirr.d data for plotting and calls IPIOG. 

IPLOG: This forms a suoroutine of OUTPT. OUTPT calls 

subroutine IPLOG to make plots of the data requested. After 
plotting, control is returned to OUTPT again for further 
decision. 

Sampling Pate ; 

For getting the sampling rate of 250 samples /sec. 
the signal has to be sampled at 40 msec, interval. The 
procedure for achieving this sampling rate is described 
below. 

While reading an analog input, the flow chart shown 
in Fig. 3.10 has been used. The time- elapsed between two 
sample is equal to the time taken by the program to 
execute the instructions from A to 3 and back to A. This 
time is found to be 39.85 msec, by trial and error. To 
obtain the required interval of 40 msec, some redundant, 
instructions have to be introduced in the return path. 

These instructions are 

LOOP MDX 2 -1 

BSO LOOP 
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Tii6S6 "two ins L.i'"u.G't ions "to-ks 10 j.issc, foi* osscu.'fcion. Tlisss 
instructions are executed repeatedl- to oltain tlie desired 
delay of .15 msec. 

I’ew s3-lient features of the data logging program are; 

(1) Kumher of points to be recorded on the magnetic tape 
is variable and can be entered through the data switches. 

(2) Program can be kept in standby position for any 
amount of time, during which experimental conditions at 
the transmitting end can be changed. This avoids the 
starting of the system again and again during data logging. 

(5) There always exisis the possibility that the signal 
amplitude ma37 incidently exceed the limit of /JDO (+ 5 Volts) 
and in this case the ADC will not show the appropriate value 
of the signal. Before recording any data on the tape 
DTIGI checks for the limit (+ 5 volts). If the signal 
amplitude is more than + 5 volts, the data is treated to 
bad and not put on the tape. Thus DTIGI recognises the 
bad and good sample. 

(4) Any record on the magnetic to,pe which is suspected to 
be bad or gives trouble while reading can be erased and 
fresh data can be logged thereafter. 

Program for Analysis of the Data ; 


Initially it is thought that the analysis can be done 
by IBM 1800, and a program is developed for ..processing the 
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data by IBM 1300. This program has been given in Appendix B. 
But later it is found that IBM 1800 gives dimensionality 
problem even after making many links and using small 
subroutines and takes more time as compared to IBM 7044. 
Program given in Appendix B works well for limited number 
of auto -correlation points. Since IBM 7044 is much more 
faster than IBM 1800, and has large memory, a program ’PROGS* 
for analysing the data by 7044 is developed. This is given 
in Appendix G. The program ’PROGS' has eight subroutines and 
few system subroutines. 

Main program PROGS READTP 

TRUAB 

AUTO 

PSD 

SMUTil 

GRAPl 

I3I'31L 

RPITGE 

PROGS s Main program PROGS dictates the input/output scheme. 
PROGS can bo operated in two modes. Pirst mode is the 
Sontinuous mode, in the sense tha.t every data on the 
input magnetic tape is analysed, and the second mode is 
the Selective mode where only selected data are analysed. 
PROGS can give paper output, punch output or output on 
magnetic tape according to the user’s requirement. In this 
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case the out nut is mainly taken on magnetic tape which is 
later brought back to IBM 1800 for plotting the result. 

RBABTP ; Since the word format of IBM 1800 and IBM 7044 
are different, conversion is needed for reading the IBM 1800 
magnetic tapo for IBM 7044. Subroutine EEADTP establishes 
the linkage between IBM 1800 and IP.M 7044 and thus EEADTP 
is very essential subroutine without which processing of 
IBM 1800 magnetic tape on which data was logged is not 
possible . 

T READ ; Subroutine TREAD brings required number of words 
from the magnetic tape to the computer memory. TREAD does 
so by cross referencing the subroutine EEADTP. 

AUTO ; Subroutine AUTO uses simple algorithm (Eq. 2.6) 
for computing the auto -correlation of the data points. 

PSD : This subroutine calculates spectral d ensity of the 

signal by taking the Eourier Transformation of the Auto- 
correlation function. Eqn. (2.15) has been used in writing 
the computer program. 

SMUTH : This computes the smooth power spectral density 

of the signal. Eqn. (2.17) has boen used in writing the 
computer program. 

ORAEl : This subroutine is used for plotting the output 

Discrete plotting of points drawn by IBM 7044 is not as good 
as the continuous plot drawn by IBM 1800. A plot is taken 
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©n IBM 7044 to just check the output, 

ISEEL ; Before recording the final output on magnetic tape 
all values are converted to integer foim and scaled such 
that no values exceeds 3 2J76® (limitation set by IBM 1800). 

This infeger transformation is re(^uired because fhe 
program developed to read this tape in IBM 1800 handles 
only integer data on the tape ancl moreover, integer mathe- 
matics is faster and simpler to deal with. 

; Output in the form of punch card can be had by calling 
this subroutine. 


Main program PROGS is very flexible and gives lot of 
options for the user. During the course of computation it 
is possible to know that upto which serial number the 
computation has been done. 

Plotting the IBM 7044 Output by IBM 1800 Plotter ; 

The program for plotting IBM 7044 output by IBM 1800 
has 3 subroutines. Program is given in Appendix D. 

Main Program P7044! RITAP 

TREAD 

PLOT 


TSA7E 
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P7Q44 ‘ Foi* one experiiaeiitial conditiioii 4 op 5 records are 
selected ano "blioir eu’to—correla’tioiis and smoo'th. power 
spectral density are averaged over 4 or 5 records and 
the averaged value of the auto-correlation and smooth 
power spectral density is plotted. This is done to compare 
the average value and the value olotained by taking singl.b 
record. 

RITAP ; It has been mentioned earlier that the word format 
of IBM 1800 and IBM 7044 is different, once again a link 
has to be established in IBM 1800 and IBM 7044, Subroutine 
RITAP serves this purpose, 

TREAB ; This subroutine brings the required number of words 
from magnetic tape to the memory of IBM 1800 by cross- 
referencing the subroutine RITAP. ' 

PLOT : The subroutine PLOT when called by P7044 draws the 
required plot and returns the control to P7044. 

TSAVB ; At the end of run, subroutine TSA7E \niloads the tape 
and requests the user to save it, 

5,4 Bxnorimontal Procedure i 

Experiment is performed in two pa.rts (a) with gas 
bubbles, (b) with nucleate boiling bubbles. 
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( a ) G-as Bubbles .^ 

In tliis casG the bubhles are gencrp.ted by passing 
air through, a bene copper tube with a hole of 1.5 dia 
in it. The copper tube is dipped in the glass jar having 
demineralized wat-iir. Rate of bubble fonua,tion is controlled 
by adjusting the air flow ra.te through the copper tube. The 
whole set up is alligned such that the collimated beam of 
light, after intercepting the bubbles forming in the glass 
ja.r, falls on the photo multi pliejc tube. The PMT senses the 
fluctuation in the light intensity and converts it to the 
electrical signal. After proper amplification the signal 
is coupled on the transmission medium and at the receiver 
end, the signal is fed to the input interface and filter 
where unwanted frequencies arc cut off. Output of the filter 
is connected to the ADC of IBM 1800. Bata acquisition is 
initiated after appropriate selection of the record length 
and sampling frequency. The digitized data are recorded on 
the magnetic tape. The few data are processed on IBM 1800 
and the majority of the data are processed on IBM 7044. 

Auto -correlation and smooth power spectral density are plotted 
on IBM 1800 plotter. 

Ohoice of Parameters; 


Bubble multiplicity and frequency of emission is 
controlled by varying the gas flow Jirate through the copper 
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tulDe. Tlie ©ffec b oi surfece ‘tension is studied "by adding 
soap solution to the water. However, only the qualitative 
study of the variables are done in "the iDrosent work. 

( b ) Nucleate Boiling ; 

In the nucleate boiling experiment, the bubbles 
are formed in the water tank on the tip of the copper rod 
heater. Rest of the procedure is same as the gas 
bi^bble experiment . 

Choice of Parameters ; 

1) Degree of Sub-oooling: 

Variation in the degree of sub-cooling is achieved 
by varying the bulk temperature of the water in the water 
tank, and the bulk temperature variation is obtained by 
vaiying tho vo It o-ge to the auxiliary heaters. 

The voltage applied to the auxiliary heaters is 
incremented in steps of 40 ■volts, the lower limit being 
140 volts and the upper limit being 220 volts. Thus four 
different bulk temperature is obtained. 

ii) Heat Rlux to the Conner Rod Heater ; 

For each setting of the bath temperature four 
different heat flux conditions are obtained for the copper 
rod heater by varying the voltage using r. variac. Minimum 
applied voltage to the copper rod is 50 volts and maximum 
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applied voltaS® is 90 volts. The incroment in voltage is 
done in steps of 20 volts. 

ii i ) Surf ac e leps io n : 

To study tho effect of surface tension on the auto- 
correlation power spectral density, soap solution is 
added to the water (Bulk temp; 85. 5^0, voltage applied to 
copper rod heater 90V). However, only qualitative 
study is done in the present work. 



CHAPTER 4 


EBSUITS 

4.1 Pat a Reduction and Presentation of Results ; 

THe correlation technique has heen used i?o determine 
the auto-correlati on and power spectral density of gas huthles ■ 
and nucleate boiling bubbles. The auto -cor relation has been 
calculated at 1000 points. The auto-correlotion is normalized 
to zero delay time. It was found that by taking 4 or 5 records 
for a,veraging tho effect of noise has been diminished and the 
accuracy is improved considerably, formalized auto-correlation 
and smoothed power .spec tral density are averaged over 4 to 5 
records and the average value of the functions are plotted • 
against delay time and frequency respectively . 

4.2 G-as Bubble ; 

Pig. 4.1a and 4.1b shows the auto-correlation and 
smoothed power spectral density plot for single gas bubble. 

denotes the normal surface tension of the water. Pv;aks 
occur at the frequency of generation, of the bubbles, and at the 
harmonics of that fundamental frequency. In Pig. 4.1b the 
first peak occurs at 5.1 cps. The amplitude of the second 
hanaonic is very small compared to the fundamental frequency. 
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Figs. 4.2j 4.3 and 4.4 show ihs C'ff'Gcfcs of surface 
tension variation on auto-correlation and spectral density. 

Here the surfaces tension o’ 2 _» <^2 increasing 

order of magnitude. It is observed that as the surface tension 
is increased moro and more peaks appear in spectral density 
curve towards the hi^r frequency. 

Many records for single bubbles are analyzed and it is 
found that when the flow rate is small, single peaJks charac- 
terizing only one distribution of bubble size appeared in 
the spectral density curve. Figs. 4. 5a and 4.5b show the 
auto-correlation and smoothed power spectral density curve 
for multibubbles when the air flow rate through the copper 
tube is increased. Here the peaks occur cl different frequencies 
with considerable amplitude indicating the splitting of tho 
bubbles with increased flow rates. 

The total amea under each of the spectral density 

curve is physically related to the power content of the signal. 

Many spectral density curve with various flow rate are 

curve 

analyzed and the area under the spectral density/arc calculated. 
It is found that area under the spectral density curve 
increases as the flow rate increases. Unfortunately, 
quantitative results c f^uld not be given since the flow rate 
measurements are not done. 
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4 . 5 Nucleate Boiliriig :- 

Pig. 4.6 to 4.8 show the auto-correlation and spectral 
density plots for different degree of suh-cooling and diffe- 
rent heat fluz to the copper rod heater. Surface tension 
effect ha-S been shown in Pig. 4.9. 

In the nucleate boiling case, the size of the 
bubbles are too small for individual detection and the 
main observation is due to the convection currents. Since 
the size is small one expects spectrum to show peaks at bi.gh 
frequencies but however^the signal is wealc compared to the 
convective currents signal which is of low frequency. The 
peaks at low frequencies are due to the convective currents. 
Analysis becomes difficult at these low frequencies and one 
has to usG frequency modulation to record such signals. 



CHAPTER 5 


OOHCLUSIOFS AHP SUG&ESTIOIS FOR FUTURE WBK 

5.1 Conclus ions ; 

The correlation technique has yielded auto-correlation 
and powor spectral density of the gas bubble in reasonable 
s-greement with the expected results. As the flow rate is 
increased bubble frequency increases and the area under the 
spectral densitj!" curve also increases '.howing the more 
power content associated with the signal. As the surface 
tension is increased more and more peaks appeared in the 
spectral density curve towards the higher frequency. 

In the case of nucleate boiling, the size of the 
bubbles emerging from the copper rod heater are too small 
and so the signal amplitude is very small as compared to 
the convective currents. In this case the peaks occurred 
at low frequency due to convpctive currents and th© analysis 
became difficult. In the nucleate boiling, it was not 
possible to get the single bubble from the copper rod heater, 
and hence the behavior of the single buoble could not be 


studied . 
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5.2 Sugg e stions for Future Work : 

1* Generation of" single bubble is a problem in the 

case of nucleate boiling, and a better method for generating 
single bubble should be used. 

2. better method for detecting the single bubble 
should be used, because the bubbles generated in nucleate 
boiling are too small. 

3. Filaments where several sites occur detection 
should be selectively done for one site. 

4. In presence of the strong convective currents, 
where the signal amplitude is comparatively low, frequency 
modulation method for recording the signal should be used. 

5. Adjustable slits should be placed on both sides 
of the tank in order to collimate the light beam at the 
desired width. 

Possible Application of the Findings ; 

Detection of the incipience of nucleate boiling 
in the reactor core can be done remotely using this tech- 
nique. This could be accomplished by placing the light 
source at one end of the coolant channel and PMT or some 
other detection probe at the other end. 










OOOTROI CARDS ROE ESECUTIOIT OR PEOGS - IE I3M 7044 


Program PROGS has been put on the tape in EBIOAP 
form. Rolloiring control cards are required for executing 
the programs 

Magnetic tape U 437 is input tape 

Magnetic tape U 440 is reload tape 

Magnetic tape U 443 is output tape 

1. • Card 1 

$ JOB MS 029 7 , TBffiOS 5 , RAOB 304 5 , EBIB 

2. Card 2 

Sequence Card 

3. Card 3 

IPAUSB MT.U440-B4,U437-C7 BOTH WITHOUT RIEO,AED 
U443-0E B9 WITH RIHG, THAi'IKS. 

4. Card 4 

16 

$ IB JOB HOSOUECB 

5 . Card 5 

$RBI.O AD U04 , MS=PiTSHOE , SRCH 

6 . Card 6 

1 -5 6 -10 11-15 16-20 21-25 

LRBIN liBOT lUNlD LAST IWAIT 



LiJBIlf i Li:3!-.l of input tape, 437 

LASOT ! Latel of output tape, + 443 only for first time, 

~ 443 for rest of the run. 

= 0 OUTPl on tape is not required. 

IITI'TLD ;= 1, unload tta.e tape 

= C, unloading not required 
LAST ; Last record on the output tape if 
LAST = 0, it is taken as 52000. 

I'-JAIT : 0 

7. Caxd 7 

1 - 5 6 -10 11-15 20-50 51-55 

JSSR HLAST I AT VI. JSIf M 

COl-TTIMUOUS MODE ; For continuous node of operation two 
blank cards are needed after Card 6. 

SELECTIVE MODE ; Card 7. 

JSBPl ; Serial number of the data to be processed. 

ILAST i Serial number of the last data to be processed. 
lETVL s Interval by which the serial number of the 
data has to be incremented. 

JS¥ ; Switches for selecting the output. 

Type of Output for Signal; 

1 in column 20 for graph 
1 in column 22 for paper output 
1 in column 24 for punch output 



CONTROL C/jRI)3 REQUIRED EOR EXECUTION OP P7044 
I3II 1800 


1. Card 1 

// JOB 

2. Card 2 

// XBQ P7044 PX 

3. Card 3 

1-5 6-20 

L7ESL IDATE 

LiJ3EL : label of the tape in 15 format 
IDATE : Month, Date, Year in 315 format 


e.g. 11 

29 

Card 4,5 


1-20 

21-40 

L/EZ 

LABY 


LZBZ : X-axis label and auto-correlation/spectral density 
LABY ; Y-axis label and auto-correlation/spectral density 

5. Card 6 

1-5 6-10 11-15 

ISTRT NLAST lEWTl 

ISTRT ; Storting serial number 

NIAST : last serial numper upto which averaging has to 


be done 



1C3¥TI : 1, read title card for auto-corrolation and 

spectral donsit 3 ^ in 30A2 format in two cards 
= 0, no title card 

6. Card h 

Same as card 6 

if ISTRT = C, program calls CIII EXIT. 






